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Abstract 

00 \ We examine the predictions for both the spin dependent and spin independent 

•O direct detection rates in a variety of new particle physics models with dark matter 

^ ■ candidates. We show that a determination of both spin independent and spin de- 

l } , pendent amplitudes on protons and neutrons can in principle discriminate different 

candidates of dark matter up to a few ambiguities. We emphasize the importance 
of making measurements with different spin dependent sensitive detector materials 
and the need for significant improvement of the detector sensitivities. Scenarii where 
' exchange of new coloured particles contributes significantly to the elastic scatter- 

^ ' ing cross sections are often the most difficult to identify, the LHC should give an 

^"'i indication whether such scenarii are relevant for direct detection. 

OT 

i — '! 1 Introduction 

■ Unraveling the properties of a new stable cold dark matter (CDM) particle is a challenge 
for ongoing or future astroparticle and collider experiments. The most convincing evidence 
£D '. for CDM so far is provided by WMAP [T] and SDSS [2J. Their precise determination of 
the relic density of CDM strongly constrains the parameter space of the various new 
particle physics models(NP) [31 HI El El El El El HTJ] - This single observable, VtcDuh 2 , 
is however not sufficient to pin down the properties of CDM even when assuming that 
O ■ the candidate is a weakly interacting massive particle, x- Additional information on 
>• ! the nature of dark matter could also be obtained from measurements of detection rate in 
different detectors, observations of a signal in photons, antiprotons, positrons or neutrinos 
produced after annihilation of dark matter and discovery and measurements of properties 
of new particles at colliders. 

Several models for new physics containing a CDM candidate have been proposed in 
the past [llj. The most popular examples of new stable weakly interacting particles at the 
electroweak scale include the neutralino in supersymmetric models [T2J [T3] , right-handed 
neutrinos pU [15j [16], scalars or vector bosons in extra dimension models p21 dU [19] . 
vector bosons in little higgs models [201 El] and scalars in extensions of the SM [221 El] • 
Right-handed sneutrinos as a CDM candidate have also been revived lately [251 126] . 
Predictions for signals in direct PH EH EE! M> M> EI] , indirect [321 E31 El ESI ESI EH 
EH1 EH1 HOI IH1 1121 1131 HI] or collider experiments have been made within each of these 
models [1H1 SSI d SZl SHI SHJ EQl EU E21 ESI EH EHl ESI [57]. Furthermore in specific case 
studies, in particular within supersymmetric models, the prospects of determining the 
properties of the new particles and from there infer a "collider" prediction for the relic 
density or for the detection rates were analyzed [SHI EH1 SSI ED]- While colliders and in 
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particular the LHC have a good potential for discovering and identifying new particles 
present in various extensions of the standard model, direct detection experiments (DD) 
will be the ones to provide evidence for a stable relic particle [611 162] . Furthermore in 
some cases, direct detection experiments have better discovery prospects than the LHC. 
The best known example is the so-called focus point region in constrained supersymmetric 
models [nSHHHESj- We therefore concentrate here on direct detection aspects and consider 
only models which offer the best detection prospects, those with a weakly interacting 
particle at the electroweak scale. 

A number of experiments are currently searching for CDM by measuring the elastic 
scattering rate on nuclei in large detectors. Their sensitivity is being improved and upper 
limits are updated regularly. The best upper limit on the proton- % spin-independent (SI) 
cross section has been recently obtained by Xenon, < 4.5 x 10~ 8 pb for a CDM of 
30GeV [66] and CDMS, < 4.6 x 10" 8 pb for a CDM of 60GeV [SI]. These limits 
are already putting constraints on the parameter space of new physics models. Limits 
on spin-dependent (SD) cross-sections are much less restrictive. The best limits are now 
obtained by KIMS for protons, <r~P < 1.6 x 10 _1 pb [HE] and by Xenon for neutrons, cx^f < 
6 x 10~ 3 pb [SD] . Indirect detection of neutrinos coming from CDM annihilation in the Sun 
sets a limit on cr^P, the best limit is from Super-Kamiokande, < 4. x lCT 3 pb [70] . 
These do not yet allow to test the most popular NP models 0. 

One difficulty in extracting precise information from an elastic scattering rate on 
nuclei is that the rate depends not only on the details of the particle physics model but 
also large theoretical uncertainties are introduced by the CDM velocity distribution, the 
nuclear form factors and the coefficients that describe the quark content in the nucleon. 
The former can be eliminated by taking ratios of rates in different materials while a 
large part of the uncertainty from the quark content in the nucleon will drop out when 
taking ratios of proton to neutron amplitudes. Fortunately many of the detectors set 
up or planned use different materials thus can be sensitive to different combinations of 
proton and neutron amplitudes. The procedure for extracting in a model independent way 
the amplitudes for spin dependent interactions on protons and neutrons was discussed in 
[73j [74]. The spin independent interactions on the other hand are basically sensitive to one 
combination of neutron and proton amplitudes, this is because all heavy materials have 
roughly the same ratio of protons to neutrons [75] . Our goal is to see what would be the 
prospects for determining the properties of CDM particles after a signal has been observed. 
We will use directly the event rates or assume that the spin (in) dependent proton and 
neutron amplitudes have been extracted. We assume that future spin independent and 
spin dependent detectors will have sufficient sensitivity to measure a signal. We will 
consider the maximal achievable sensitivity to be ~ 10 -10 pb, for example with 
Eureca [TB] and cr^f « 4 x 10~ 7 pb, as in the COUPP proposal [77] . Although this last 
value requires a significant improvement in SD detectors, we emphasize the importance 
of the SD interactions in determining the properties of the CDM candidate. 

Comparative studies of the prospects for direct detection in new physics models have 
been performed in [TSJ [75] . Recently a comparison of the SI detection rates and rates 
for indirect detection of neutrinos in the case of the MUED, little Higgs and MSSM 

1 Note that DAMA/LIBRA have very recently confirmed their annual modulation signal [7T|. We will 
not consider this result as it seems to be incompatible with other searches unless the CDM particle is 
below lOGeV [72], in the sample models we consider CDM candidates are rather in the 30-1000 GeV 
range. 
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models was presented [HO]- The potential of a combined measurements of SI and SD rates 
to distinguish MUED from MSSM with COUPP using two different materials was also 
examined [HI]. We expand on these analysis in many ways. First we examine a larger 
class of MSSM models, second we rely heavily on detectors sensitive to SD interactions 
and third we insist on the importance of using different materials to extract both the 
neutron and proton amplitudes. We also take into account uncertainties from the quark 
coefficients in nucleons and use an improved calculation of the direct detection rate [83] . 
This is a first step towards a more general analysis where one would combine information 
from both direct and indirect detection as well as from collider searches, see also [79] . 

This paper is organized as follows: after setting up our notation in section 2, we 
summarize in section 3 the predictions for the SI and SD cross sections on nucleons as 
well as for the ratios of SD and SI amplitudes on protons and neutrons in different CDM 
models. We consider Majorana fermions (in particular the neutralino in SUSY), a right 
handed Dirac neutrino, as well as vector and scalar particles. The results of our scans 
over the parameter space for each sample models are presented in Section 4 where we 
show which models can in principle be distinguished by measurements of both SI and SD 
amplitudes on protons and neutrons. The predictions for the scattering rates on various 
nuclei are then compared in section 5. Finally in section 6 we briefly mention the case 
where a signal can be observed only in the SI interaction. Our results are summarized in 
section 7. 



2 Direct detection 

The total scattering cross section of a DM particle, \i °ff a point-like nucleus for spin 
independent interactions reads 

^^ rZ + UA-Z)f (1) 

7T 

where \i x = m x MA/{m x + Ma) is the reduced x- nuc l eus niass and Ma the mass of 
the nucleus. Te proton(neutron) amplitudes, are related via some coefficients to the 
amplitudes for %-quark scattering, X q . For example for scalar interactions of Majorana 
fermions, in the notation of [H3] 

K,n = £ F q ^ , (2) 
9=1,6 

where f^ ,n describes the contribution of quark q to the mass of the nucleon. The quark 
coefficients for scalar interactions have large uncertainties [84J. To take these into account 
we vary the input parameters of micrOMEGAs2.2 [H3] in the range 

a nN = 55 - 73MeV and a = 35 ± 5 MeV (3) 

which in essence amounts to varying the s-quark content in the nucleon in the range 
0.19 < fs < 0.56. The heavy quarks coefficients, fg, are related to those of the light 
quarks [H5]. I n the case of a Dirac fermion with an effective vectorial interaction, the 
coefficients that describe the quark content in the nucleon just count the number of 
valence quarks and therefore have no theoretical uncertainty [85] . 
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For spin dependent interactions, the point-like nucleus cross-section reads 

where J a is the total spin of the nucleus and S^„ are obtained from nuclear calculations. 
The SD nucleon amplitudes, £ Pjn □ are related to the quark amplitudes, 

= A( i p,n ^ ( 5 ) 

q=u,d,s 

where the coefficients Aq p,n have been estimated for light quarks [85] 

A£ = 0.842 ±0.012; A^ = -0.427 ± 0.013; Af = -0.085 ± 0.018 

K = *Z\ K = K\ a: = a? (6) 

In the numerical analysis we will allow the coefficients to vary within their la range. 

The recoil energy distribution measured in a detector further contains some depen- 
dence on the nuclear form factors as well as on the CDM velocity distribution. 

dN 2M det t po 



dE 7r m x 



[Fl{q) (X P Z + X n (A - Z)f (7) 

+ i^j^{sMa + s m {q)Ui + SiMiT) HE) (8) 

where £i = £ p ±£ n and £o = £p — Cn- Fa(q) is the nuclear form factor for scalar interactions 
and Sij(q) the form factor for spin dependent interactions, both depend on the momentum 
transfer, q = [82]. p is the local neutralino density, M^t the detector mass and 

1(E) the integral over the velocity distribution 

1(E) = f°° ^dv (9) 



mi n 



where v min = (^^-) 1 ^ 2 - To compute the cross sections on nucleons and the event rates 
we rely on micrOMEGAs2.2 [83]. 

All information on the CDM model is contained in the amplitudes X p , A n , £ p , ^ n as 
well as in the mass of the CDM. Once a signal has been observed, one could use data 
from different detector materials to extract information on these amplitudes [86]. This 
evidently necessitates making some assumption about both the halo velocity distribution, 
the dark matter distribution as well as on the nucleon and nuclear form factors. The 
dependence on the velocity and on the dark matter distribution however drops out when 
taking ratios of the nucleon amplitudes. Furthermore, some of the uncertainty from the 
quark coefficients in the nucleon also drop out. This is because the sea quarks coefficients 
which give the dominant contribution are identical for protons and neutrons. 

We choose as independent parameters the ratios (f>, £/A+, X p /X n and A + which charac- 
terizes the overall SI rate, 

£ p = £sin0; £ n = £cos0; A+ = A p + 1.4A n . (10) 



J Note that our definition of the nucleon amplitudes differ from the usual convention where one uses 
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Note that the factor of 1.4 in A + depends on the ratio of protons to neutrons in the nucleus. 
Typically this ratio does not vary much [75], our choice gives the maximal sensitivity in 
heavy nuclei. Nuclei that are sensitive to SD interactions do so primarily through an 
unpaired nucleon, this means they have either S p or S n 7^ and have little sensitivity to 
the interference term £ p £ n in eq. H] [H7J. The sign of = atan(^ p /^ n ) is therefore hard to 
determine. 

The mass of the CDM candidate can be determined from the nuclei recoil energies. 
This works best when M x « 100 GeV [88J although a new method to improve the mass 
determination for a heavy DM particle using signals from two different detectors was 
proposed recently [89] . 

3 DM models 

We consider a selection of models representative of different CDM candidates in the 
30GeV-lTeV range: Majorana fermion (the neutralino in the MSSM), Dirac fermion (a 
right-handed neutrino), gauge boson (the heavy photon in little Higgs models or the B\ 
in MUED models) or scalar particles in extended Higgs models. The predictions for the 
rates for direct detection have been studied in all these models and rates can vary by 
orders of magnitude within each model [7J [TT]. For each CDM candidate the dominant 
process for elastic scattering influences the overall scattering cross-section as well as the 
relations among the proton/neutron amplitudes. We will explore these relations within 
models representative of each type of CDM. A summary of the different mechanisms for 
CDM elastic scattering in various models is provided in Table [U Note that two special 
subclasses of the MSSM have been introduced. The main difference between these two 
classes is the range of mass of squarks, in MSSMH they are heavy and therefore do not 
contribute to DD. 



Table 1: Dominant mechanism for CDM-nucleon elastic scattering 



Model 


CDM 


Nature 


SI 


SD 


MSSMH 


Neutralino 


Majorana fermion 


Higgs 


Z 


MSSMQ 


Neutralino 


Majorana fermion 


Higgs+squark 


Z + squark 


RHN 


Vr 


Dirac fermion 


Z +Higgs 


Z 


MUED 


B x 


Vector boson 


Higgs +KK-quarks 


KK-quarks 


LHM 


A H 


Vector boson 


Higgs + Quarks 


Quarks 


IDM 


H" 


Scalar 


Z+Higgs 





3.1 MSSM 

In the MSSM the CDM is a Majorana fermion, the neutralino, Xi The nature of the 
neutralino, whether it is mostly bino or contains a mixture of Higgsino or wino, strongly 

3 We do not consider the case of the sneutrino CDM which usually gives too large detection rate unless 
its coupling to the Z is suppressed. [25l ESI [90] 
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influences the annihilation mechanisms and the CDM relic density. For direct detection 
one gets two types of contributions, Higgs and squark exchange for SI interactions and Z 
and squark exchange for SD interactions. In general the Higgs and Z exchanges dominate 
since the squark contributions suffer from a mass suppression (the squarks are generally 
at the TeV scale). We will consider two categories of MSSM models. In the first, MSSMH, 
sfermions are heavy (2TeV) and do not contribute to DD. In this model CDM annihilation 
requires a lightest neutralino with some Higgsino or wino component. In the second, 
MSSMQ, we force one of the squark masses to be Mg L R < 2M X . We introduce these 
classes of models as a way to quantify the impact of the squarks in DD. Note that these 
two types of MSSM models can be easily distinguished at LHC which can probe the 
squark sector up to more than 2TeV [HH E2] • In both cases we will assume the sleptons to 
be heavy since sleptons do not contribute to direct detection. However one should keep 
in mind that sleptons can contribute to the CDM relic density, both through annihilation 
or coannihilation, so our analysis is not completely general when confined to models that 
are in agreement with the measured value for the CDM relic density. Such light sleptons 
can be searched for at the LHC. 

Taking into account the dominant Higgs exchange diagram only, the spin independent 
interaction reads 

A^v = -™ NJ J^- £ + if + f t N )9k lU u + (f d N + f? + f b N )9k ldd ] 9h iXX \ (11) 
AM w c w m h . 

where hi = h,H, g huu = cos a/ sin/3, g hdd = - sin a/ cos/5, g Huu = sin a/ sin/3, g Hdd = 
cos a I cos/3 and a is the Higgs mixing angle. In the decoupling limit, at large M#, 
sin a = — cos (3. The SUSY-QCD corrections can shift the Higgs couplings to down-type 
quarks, especially at large values of tan (3. These corrections are taken into account in 
the numerical analysis bur for simplicity will be omitted from the discussion here. The 
couplings of the light Higgs to Xi reads 

9h X x = (cosaZ 14 + sinaZ 13 )(c w Z 12 - s w Zu) (12) 

where Z±j describe the field content of the LSP [92]. Clearly for the Higgs exchange to 
contribute requires a LSP with some Higgsino component (Z 13 , Z 14 ^ 0). In all models 
where Higgs exchange dominates we expect X p = X n within a 2% accuracy. This is because 
the the quark coefficients in protons and neutrons are the same for heavy quarks and the 
largest coefficient is the one for s-quark. 

For the spin dependent amplitude, the Z exchange contribution reads, 

e w = -i ( A f -Ay-Af) Ij ^(^-^) (i3) 

When the squark contribution is negligible the ratio of proton to neutron amplitudes 
£p/£n is therefore totally independent of the neutralino coupling to the Z. We expect 
tan0 = £ p /£ n ~ —1.14 ± 0.03 when considering the range for the quark coefficients 
specified in eq. El This value for tan0 is expected in any model where Z exchange 
dominates the spin dependent interaction. 

These simple relations are spoiled in models where squarks are light unless squarks of 
different flavours are nearly degenerate in which case we still expect \/\ n ~ 1. Strong 
corrections to this ratio can be found when m,q L R m m x since in this case twist-2 operators 
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give a large contribution that can even cancel the leading squark contribution [S3]- Note 
that a twist- 2 contribution, being proportional to the quark hypercharge, is larger for 
u-squarks so will contribute mainly to the proton amplitude. In principle the impact 
of light squarks is more important for spin dependent interactions because of a possible 
cancellation between u and d-type quark coefficients. However the squark exchange is 
dominant for SD cross sections that are usually too small to be measured even in ton- 
scale detectors. 

The ratio of SD to SI amplitudes also characterizes the model. When sfermions are 
heavy, the relic density of dark matter favours a LSP with some wino or higgsino compo- 
nent. In such models the ratio of spin dependent to spin independent interactions depends 
strongly on the higgsino component of the LSP, eq. [TT] and eq. H31 and predictions can 
vary in a wide range. 



3.2 Right-handed neutrino model (RHNM) 

A model with warped extra dimensions where the CDM is a right-handed Dirac neutrino 
was proposed by Agashe and Servant [93] . This model contains both new fermions and 
new gauge bosons at the (multi-)TeV scale which interact mainly with third generation 
fermions. This model can be used as a prototype of a more general class of models 
with a right-handed Dirac neutrino as CDM [93]. Whether or not there are additional 
quarks or gauge bosons, because of the large mass scale involved, the most important 
contribution to elastic scattering of the right-handed neutrino on nucleons is due to Z and 
Higgs exchange [93J. What is peculiar in this class of models is that the CDM is not a 
Majorana particle, so there is an important contribution of Z exchange to both SI and 
SD nucleon scattering [83]. Typically in this model the elastic scattering cross sections 
are large and direct detection poses one of the strongest phenomenological constraint on 
the model [93] . 

For the dominant Z exchange contribution to the spin independent interactions 

_ g?e(l - 4^) _ gje 

8M 2 z s w c w a An 8Mls w c w [ } 

where g^ R is the parameter that describes the coupling of vr to the Z. This coupling is 
induced through mixing so is suppressed with respect to the SM couplings. The neutron 
and proton amplitudes are directly related 



^ = (1-44)^0.09 (15) 

for Syy = 0.228. If the vr also couples to the Higgs, both \ p and A n will receive the same 
additional contribution, thus modifying the simple relation, eq. HS1 In the numerical 
analysis we will include a generic coupling of the Higgs to the neutrino gjj [94] . 
For spin dependent interactions which also proceed through Z exchange, we get 



b'\* where W= ^LJl1 



6r = 2J b>A» where b> = " **L_ ^ 

q=u,d,e ' 



and g R — gf = —(g R — g^) = 2s e c ■ As for the MSSM, the ratio of proton to neutron 
amplitudes, tan</> = £ p /£ n = —1.14 ± 0.03. 
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In the limit that the Higgs contribution is negligible, the ratio of SD to SI amplitudes 
is also independent of the details of the model with £/A + = 1.06 ± 0.02 when varying the 
quark coefficients in the range specified in eq. El 

3.3 Universal extra-dimensions model (MUED) 

In the universal extra dimension model (UED) potential dark matter candidates include 
a KK gauge boson, a KK neutrino, a KK scalar or a KK graviton (95j [TQl [36] . We restrict 
our analysis to the minimal UED model (MUED), in which case the CDM is either the 
first KK level of the hypercharge gauge boson, B 1 or the KK graviton. We will consider 
only the former possibility since the graviton has small direct detection rates. CDM 
scattering on nucleon proceeds both through Higgs exchange and KK-quark exchange. 
For spin independent interactions, the nucleon amplitude reads [T9] 



where the sum is over all quark flavours and g x = e/cw, Y qL = 1/6, Y UR = 4/6, Y UR = 
—2/6. The first term arises from Higgs interactions and the second term from KK quarks 
exchange. In this model it is quite natural to have a large contribution from KK-quarks 
since they are nearly degenerate with the CDM. We include radiative corrections to level 
1 KK states [18] which lead to mass splittings between KK quarks and B 1 . Note that the 
Higgs contribution is suppressed compared to the MSSM by a factor m w /m Bl as well as 
by the Higgs mass which is usually larger than in the MSSM. Nevertheless one expects 
\/ A n ~ 1 as in models where the Higgs exchange dominates because all new quarks are 
nearly degenerate. 

For spin dependent interactions, the amplitude reads 



and is solely due to KK quarks exchange. One can easily show that tan^> = £ p /£ n ~ —3.5 
independently of the parameters of the model as long as all KK quarks are degenerate. 
The ratio £/A + can be large and is controlled by the B\ mass and by the mass splitting 
with the KK-quarks when these dominate the SI interaction. 



In the little Higgs model with T-parity, the dark matter candidate is the lightest new 
heavy neutral gauge boson Ah [96l [97] . This model therefore shares many aspects of the 
MUED model just discussed, the CDM is a gauge boson, spin independent interactions 
are due to Higgs and heavy quark exchange, while only the latter contributes to spin 
dependent interactions. The expressions for both SI and SD amplitudes are the same as 
above, eqs. fTTfTSl There are however two important differences between these models: 
first the hypercharges of the heavy quarks are small, Y qL = 1/10, Y qR = 0. Second the mass 
splitting between the new heavy quarks Q and Ah is typically much larger than in the 
MUED model, which means that the heavy quark contribution to DD is suppressed. One 




(17) 




(18) 



3.4 Little Higgs model (LHM) 



S 



therefore expects an overall low rate a^ N and X p / \ n ~ 1 when Higgs exchange dominates 
or heavy quarks are nearly degenerate. Because there is no Z exchange diagram the SD 
interaction should also be much suppressed unless one artificially requires a small mass 
splitting between the heavy photon and heavy quarks, eq. [HI It is only in this case that 
one expects to have a detectable cross section. Then \<j)\ will depend strongly on the mass 
difference between the heavy photon and the lightest new quark and should be large. 

3.5 Scalar dark matter 

Simple models with an additional scalar field that is basically decoupled from the SM 
sector have been proposed [221 [231 [21] • In these models, the CDM candidate is a new scalar 
field. We consider the Inert Doublet Model (IDM) [23]. a two Higgs doublets extension of 
the standard model with a symmetry. One of the two doublet and the usual standard 
model particles are even under this symmetry. The new particles of the model are a 
neutral {H°), a pseudo (A) and a charged {H + ) scalar. Depending on the parameters of 
the model, the dark matter candidate can be H° or A. Only spin independent interactions 
can occur through either H°q H°q and H°q — > Aq. The latter has to be kinematically 
forbidden, that is Ma — M^o > 100 keV, to respect the current experimental constraints. 
The h exchange cross section driven by an effective coupling A^ [23] is : 




where x = H ■ 

4 Results 

Here we present numerical results for each of our sample models. Amplitudes and cross 
sections for direct detection are computed with micrOMEGAs2.2 and in each case include 
all tree-level diagrams, the contribution of twist- 2 operators as well as QCD corrections. 
Additional SUSY-QCD corrections are included in the MSSM as discussed in [S3]- The 
computation of the CDM relic density is also based on micrOMEGAs2.2 [98j [99]. We fix 
mt = 172.6 GeV. We also restrict the parameter space to a CDM particle roughly below 
the TeV scale simply because detectors are not as sensitive to heavier CDM particles. We 
also never consider > 500 GeV, although allowed in some models such a Higgs gives 
a small contribution to direct detection. 

We first summarize for each model the predictions for both a^ N and cr^®. We always 
impose the upper limit from the relic density of dark matter Qh 2 < 0.136 |100j in our 
scans as well as other model dependent constraints on the parameters of each model. 
When specified we also impose the lower bound Qh 2 > 0.094 [lOOj . We then compare 
the ratio of amplitudes on neutrons and protons (0, £/A+, X p /X n ) before comparing the 
rates on various nuclei. We have taken into account the theoretical uncertainty in the 
coefficients that relate the amplitude for quarks to the one in nucleons by varying the 
input parameters in the range specified in eqs. I3ll6l 





(19) 
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4.1 Predictions for a^ N and orjy 
4.1.1 MSSM 

We consider two specific classes of the generic MSSM, as mentioned above. In the first 
class, MSSMH, squarks are heavy and we assume only universality among two of the 
gaugino masses at the GUT scale, that is M3 = 3M 2 at the weak scale. In the second 
class, MSSMQ, we allow for light squarks, for simplicity we also impose full universality 
of the gaugino masses, which leads at the weak scale to M 3 = 3M 2 = 6M1 . In all cases 
we assume heavy slept ons. 

For each model we have scanned over 10 5 scenarii varying the model parameters defined 
at the weak scale in the range 

100 GeV < Mi < 1000 GeV; 100 GeV < n < 2000 GeV 

100 GeV < m A < 2000 GeV; 2 < tan/5 < 52 (20) 

and 

MSSMH : 100 GeV < M 2 < 1000 GeV ; Mq L R = 2 TeV 

MSSMQ : M 2 = 2Mi; M fe(H) < 2M X , Mg R[L) = 2 TeV (21) 

For this range of parameters the mass of the DM particle does not exceed 1 TeV. In each 
case the LEP limits on Higgs and SUSY particles are imposed as well as the upper bound 
on the CDM relic density. 

The predictions for the SI and SD cross-sections in MSSMH are displayed in fig. [TK to- 
gether with the reach of future ton-scale detectors. The absolute bound for CDMS/Xenon 
is indicated only to guide the eye as this limit depends on the CDM mass. Virtually all 
scenarios will be accessible to future searches for SI interactions. This is a direct conse- 
quence of imposing the constraint from the relic density which requires a neutralino DM 
with some higgsino component for efficient annihilation. This then automatically leads 
to a Xi coupling to the light Higgs hence to a non negligible cross section for SI elastic 
scattering. The smallest cross-sections in fig. [T^l correspond to a neutralino with a very 
small higgsino fraction that nevertheless annihilate efficiently because it does so near a 
heavy Higgs resonance. Models with a relic density within the WMAP range (rather than 
only below the upper bound) almost span the full range of predictions for SI and SD 
cross sections, although many of the scenarios with the largest a^ N have a large higgsino 
component and are associated with a small value for the relic density due to the effi- 
cient annihilation into W pairs. The higgsino component also induces a coupling to the 
Z hence leads to SD interactions. Because these interactions are not coherent in several 
cases the predictions can be as low as a^f ~ 10~ 9 pb, much below the expected reach 
of future detectors. Note that because Z exchange dominates SD interactions, the rate is 
directly related to the Z\X coupling which is proportional to Zf 3 — Zf 4 , eq. [131 A rate 
measurement will therefore set a limit on this coupling assuming the MSSMH, see fig. [lb. 

In MSSMQ, the range of predictions for cr?l is roughly the same as in MSSMH, see 
fig. [2k although cross sections below the reach of future SI detectors can be expected in 
a few cases. Furthermore large cross sections for SD interactions can be expected even 
when SI ones are quite low. In general this occurs in scenarios with light squarks. There 
is no explicit correlation between a^ N and the mass of the CDM, see fig. [lb, although 
models with a neutralino around 60 GeV that annihilate near a light Higgs resonance can 
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_q _B _7 _fi _S -4 — ^ -9 -4- -9 -1 

10 10 10 10 10 10 10 10 10 10 10 10 1 

cr xp SD (pb) Z 13 2 -Z 14 2 

Figure 1: a) Predictions for vs in MSSMH. In blue the scenarios that satisfy 
the WMAP upper bound and in green those that have 0.094 < Qh 2 < 0.136. For easy 
reference the present absolute lower limit from CDMS/Xenon is indicated (dash) as well 
as future limits from large scale detectors (full) b) as a function of Zf 3 — Z 2 A in 
MSSMH. 

have a small cross section. We explicitly display in fig. [3b, the CDMS exclusion limit for 
both MSSMH and MSSMQ. Many scenarios are excluded even when taking into account 
a large uncertainty (up to a factor 3 [HI]) in the exclusion limit that could arise from the 
DM distribution. 

4.1.2 RHNM 

In the right-handed neutrino model we use as free parameters the mass of the CDM, its 
coupling to the Z, g^ R , and to the Higgs, gu, as well as the mass of the Higgs. We assume 
all other particles in the model to be above 3 TeV and therefore do not play a role in 
direct detection. We perform a scan over 100000 models varying the free parameters in 
the range 

30 GeV < m UR < 1200 GeV ; 120 GeV < m h < 500 GeV 

0.001 < g v z R < 0.01 ; 0.01 < g H < 0.25 (22) 

The range of g^ H is chosen so that the upper bound on Qh 2 is easily satisfied while not 
giving too large a^ N whereas the range for the Higgs coupling gu is set so that the Higgs 
can potentially play a role in DD. Models with a Dirac right-handed neutrino often have 
an extended gauge sector. Since this is mostly relevant for the annihilation of a CDM 
particle beyond the TeV scale, we can safely ignore this sector in our analysis. 

In this model one expects rather large rates for SI interactions, and this is in fact 
the most severe constraint on the model. Furthermore a strong correlation is expected 
between the SI and SD rates as seen in fig. Indeed these rates are governed by the 
Standard Model axial and axial-vector Zqq couplings. Note that the limit extracted from 
a^ir has to be rescaled to take into account the fact that in this model X v <C X n . In 
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Figure 2: a) Predictions for vs cr^jf in MSSMQ. In black the scenarios that satisfy 
the WMAP upper bound and in green those that have 0.094 < Qh 2 < 0.136. For easy 
reference the present absolute lower limit from CDMS/Xenon is indicated (full) as well 
as future limits from large scale detectors (dash) b) a^i as a function of the neutralino 
mass in MSSMQ (black) and MSSMH(green) 

practice it means rescaling the limit by a factor 2-3 depending on the material. Even 
taking this factor into account fig. [3h shows that models that will not be excluded in the 
near future predict a low rate for c^P- The mass of the CDM allowed in this model is 
either near M z /2, M H /2 or above 500 GeV [M]. 

4.1.3 MUED 

For the computation of the direct detection rate in the MUED model we include the 
level one KK quarks as well as the lightest Higgs exchange. We ignore the level 2 Higgs 
since in elastic scattering cross sections a heavy Higgs suffers from a mass suppression. 
Furthermore the coupling of the level 2 Higgs to the B 1 is loop induced hence suppressed. 
Note however that because M H 2 « 2M X the level 2 Higgs plays a role in the computation 
of the relic density [1011 1102] . The impact of neglecting this coupling on our analysis is 
not significant as the majority of the models already has Qh 2 < 0.136. 

The free parameters of the model are 1/R, the inverse size of the extra dimension that 
determines the mass of the KK states, A, the cutoff scale and rrih the lightest Higgs mass. 
We scan over 10 5 scenarii with the three free parameters of the model in the following 
range 

300 GeV < 1/R < 1300 GeV; 3 < AR < 50; 120 GeV < m h < 500 GeV (23) 

The precision electroweak constraints set the lower bound on |103j while perturba- 
tivity and unitarity constraints set a range for AR [104] . The mass of all KK states are 
computed including one-loop corrections |105j . The radiative corrections will induce a 
small mass splitting between the level one B l boson and KK-fermions. Such splitting is 
typically 1-2% for KK-leptons and 5 — 25% for KK-quarks and strongly influences the 
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Figure 3: Predictions for vs in MSSMQ (black), MUED (red), LHM (green), 
RHNM (pink), b) function of the CDM mass, same color code as a) with in 

addition the model IDM (black) 



direct detection rate. We also insure that the CDM relic density satisfies the WMAP 
upper bound and that the charged Higgs is not the CDM. 

The SI cross sections are suppressed by the heavy B\ mass, eq. [TTJ, the larger cross 
sections are therefore expected for the lighter CDM particles, see fig. [3b- Typically, more 
than an order of magnitude improvement in detectors sensitivities is needed to probe 
the parameter space of the model and a large fraction of the models, specially those 
with a CDM at the TeV scale, will remain inaccessible to the large scale detectors. The 
main characteristic of this model is the correlation between SI and SD cross sections, 
this is because the heavy KK-quark exchange contributes to both modes. As a result SD 
interactions could be accessible in cases where rates are too low for SI interactions. This 
is in sharp contrast with the MSSMH. 



4.1.4 LHM 



The LHM with T-parity contains in addition to heavy gauge bosons, heavy T-odd fermions 
as well as a new T-even heavy top quark. We choose as free parameters the Higgs mass, 
/, k and s a . f sets the scale of the heavy gauge bosons and fermions in particular the 
heavy photon of mass 



Ma 



v/5 



1 - 



5v z 



(24) 



with v the usual vev of the Higgs. k is an additional parameter that enters the fermion 
masses, for example for a heavy down- type quark, = \f2nf '. For simplicity we assume 
an universal factor n for all heavy fermions. s a depends on the ratio of the Yukawa 
couplings of T-even and T-odd top quarks. |106j This parameter enters the top quark 
mass as well as couplings involving standard and heavy top quarks. 
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We scan over 10 5 scenarii varying randomly the free parameters in the range 

500 GeV < f < 3000 GeV; 120 GeV < m h < 500 GeV 

0.1 <s a < 0.96; 0.11 <k<1 (25) 

We impose the LEP limits on the production of heavy quarks as well as on the Higgs 
mass. 

The rates for both SI and SD cross sections are in general quite low, even below the 
scale in fig. As we have explained before this is due to the small hypercharge of the 
heavy quarks as well as to their generally large mass. Models that could lead to a signal 
in either the SI or SD channel are those where the mass splitting between heavy quarks 
and the heavy photon is between 1-10% or slightly larger if the heavy quarks are around 
100 GeV. Furthermore the heavy photon has to be rather light with M Ah < 400 GeV, 
see fig. A Higgs near the lower LEP limit also helps increasing the signal for SI 
interactions. Note that since it is the new quarks that couple to the nucleon that need to 
be light, mainly the first and second generation, the recent Tevatron limit on the heavy 
top quark [ID 7] does not play a role here. 

4.1.5 IDM 

In the IDM, the free parameters are those of the Higgs potential [23]. We choose to use 
rather the physical parameters, the masses of the CDM candidate, m#o, the light scalar, 
m/j, the pseudoscalar, m A and the charged higgs, as well as two parameters of the 
Higgs potential fi 2 and A2. Our numerical results are not very sensitive to the value of A2 
so for simplicity we fix A2 = 0.1. Other free parameters are varied in the range 

10 GeV < m H o < 1200 GeV; 115 GeV < m h < 500 GeV; 10 GeV < /x 2 < 1200 GeV (26) 

with in addition the following range for the mass differences 

5GeV < m A - m H o < 15GeV; 40 GeV < m H+ - m H o < 50 GeV (27) 

if m H o,fi 2 < 100 GeV, otherwise 

3GeV <m A - m H o < 6GeV; 5 GeV < m H+ - m H o < 10 GeV. (28) 

We impose the following constraints on the model : vacuum stability and perturbativity 
conditions on the potential parameters, LEP limit on the charged Higgs, contribution to 
the Z boson width and electroweak precision constraints [44j . 

The rates for a^ N varies over several orders of magnitude and the masses of the CDM 
particle ranges anywhere from 50 GeV to the TeV scale, fig. [3b. Note however that once 
one imposes a lower bound on Qh 2 , the ranges for the masses and the direct rates are 
severely restricted, see section 14.41 

4.2 Discriminating models : amplitudes for scattering on pro- 
tons and neutrons 

The ratios of proton to neutron amplitudes apart from being free of large theoretical 
uncertainties provide a good model discriminator. For this to be useful one has to assume 
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Figure 4: Predictions for a) = atan(^ p /^ function of f b) £/A c)X p /X n d) M x 

in models MSSMH(blue), MSSMQ(black), MUED(red), LHM(green) and RHNM (pink). 
Only models for which a^ 1 > 10. ~ 10 pb and a^P or a?B > 4. -7 pb are included. 



that these quantities can be measured, this means that in this section we keep only models 
for which a^ N > 1(T 10 pb and > 4. x 1(T 7 pb. The more challenging case with only 
a detectable SI cross section will be discussed in the subsection 14.41 

The results of the parameter scan for the five models under consideration are displayed 
in fig. Hh-c for = arctan(£ p /£ n ) vs £, £/A + and X p /\ n . The ratio of SD neutron to 
proton amplitudes, tan 0, can discriminate models where the SD interaction is dominated 
by Z exchange (MSSMH and RHNM) from those where it is dominated by (s)quark 
exchange (LHM, MSSMQ and MUED). The ratio of SD/SI amplitudes, £/A+ which can 
be much larger in the MUED or in the MSSM could provide further discrimination. The 
parameter X p / X n can in principle disentangle further some models where the SI interaction 
is dominated by H exchange or by Z exchange (RHNM), see fig. Hfc. Unfortunately in 
practice different materials are not very sensitive to this quantity. Note also that the effect 
of twist-2 operator in MSSMQ and LHM can be important and lead to large corrections 
to the expected value of X p /X n ~ 1. 

The mass determination in DD experiment from the shape of the energy spectrum 
could in some cases provide additional information to discriminate between models. In 
particular one could distinguish the LHM, which allows a CDM in the range Ma h ~ 
50 — 120 GeV from MUED which requires a heavy CDM particle and even sometimes 
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from the MSSMQ which predicts a large range for the masses fig. HJl. The LHC, with 
its potential for discovery of coloured particles, will establish whether or not coloured 
particles could play a role in direct detection. Indeed in all LHM predicting a signal in 
DD or in the MSSMQ models, the heavy quarks and squarks can be produced easily as 
they lie well below 2TeV. The heavy quarks can be just beyond the LEP exclusion bound, 
M Qh = 100 - 400 GeV in the LHM while in the MSSMQ models that have a large value 
for |0| squarks can be as heavy as Mq < 900 GeV. The heavier squarks occur when the 
neutralino has a large higgsino content. We will not pursue a detailed analysis of what 
can be measured at LHC, this is beyond the scope of this paper. We note however that 
the mass splitting is an issue as regards the LHC potential for discovering new coloured 
particles, for small mass splitting the signals for the new particles will be hard to extract 
from the background. This could be crucial for the LHM where the mass splitting between 
the CDM and the heavy quark is below 10%. In the MUED model the mass splitting is 
between 6-22% and has been shown to be sufficient for providing a signal in four leptons 
+ missing energy channel |45j . In MSSMQ the mass splitting is also typically around 
20%. To completely discriminate between the models at the LHC would however also 
require spins of the new coloured particles to be measured |108] . 

4.3 Direct detection rates on nuclei 

Having established that a combined measurement of the amplitudes for SI and SD interac- 
tions on protons and neutrons can in principle distinguish between the underlying particle 
physics models, up to a few ambiguities, we now compare various models predictions for 
quantities that are closely related to the observables. As before we only include scenarios 
that could eventually lead to a signal in a large detector in both the SI and SD mode. We 
consider a selection of nuclei that are currently used in large detectors. Those include the 
nuclei sensitive only to SI interactions such as 40 Ar 76 Ge as well as nuclei with an odd 
nucleon that are also sensitive to SD interactions on either protons, 19 _F ; 23 Na, 127 J, 133 Cs, 
or neutrons 29 Si 73 Ge, 129 Xe, 131 Xe. For each nucleus N we compute the total event rate, 
n(N) for a recoil energy above 2 keV. For heavy nuclei, say l29 Xe, the rate is correlated 
with the value of A + while for light nuclei like 19 F the correlation is spoiled by the SD 
contribution. This is illustrated in fig. [3] for both MUED and MSSMQ. The total number 
of events varies over several order of magnitude for these two models. The rates are gen- 
erally expected to be larger for heavy nuclei, especially in models with a suppressed SD 
contribution like RHNM. In this model the total number of events in n(F) varies between 
0.4-4. x 10" 3 (events/kg/day). 

To eliminate as much as possible the astrophysical and nucleon ambiguities we com- 
pare ratios of rates for scattering on different nuclei. We define the ratios Rni/n 2 = 
n(Ni) / n(N 2 ) . We first compute the ratios Rn 1 /n 2 f° r SD proton sensitive over SD neutron 
sensitive nuclei in our sample models. Such ratios are expected to feature a dependence 
on tan <fi = £ p / £ n as well as on £/ A + when a heavy nuclei is involved. 

The results of a scan over the parameter space of each of our sample models are 
displayed in fig. The comparison of Rf/xb □ and Ri/si provides a good model dis- 
criminator with in particular only a small overlap between the predictions of MSSMH 
and MUED. This is a direct consequence of the large value for tan</> and £ in MUED. 

4 For illustrative purposes we use 129 Xe in the figures, similar results are found for 131 Xe. 
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Figure 5: Total rates (events/kg/day) a) n(F) and b) n(Xe) vs A + c) n(F) and d) n(Xe) 
vs £/A + in models MSSMQ (black) and MUED(red). Predictions for rates in LHM (green) 
and RHNM (pink) are also included in c-d. Only models for which a SI > 1. x 10~ 10 pb 
and a SD > 4. x 1CT 7 pb are included. 



In this model, 127 1 is very sensitive to SD interactions since £ p is enhanced, thus Ri/si 
is determined by the SD interaction and is large. On the other hand in MSSMH, 127 1 is 
mostly sensitive to SI interactions thus Ri/si can be reduced significantly in the scenarios 
where £ is large. Recall that those are the scenarios with a large Higgsino fraction. Such 
scenarios are precisely those that lead to a value for Rp/xeO(l) and that could have been 
confused with MUED. Indeed in MUED most scenarios predict Rf/Xb > 1 because of 
an important SD amplitude. The LHM and MSSMQ scenarios that predict large values 
for tan0, those with a (s)quark that is almost degenerate with the CDM particle have 
predictions similar to MUED for both ratios. Note that in the MSSMQ when the squark 
contribution is important there can be a partial cancellation between the various quarks 
contributions in the neutron amplitudes such that £ n <C £ p (0 ~ —n/2). Then n(F) and 
to a lesser extent n(I) are enhanced but not n(Si), these scenarios correspond to the few 
points in fig. [6^ with large R(F/Xe) and large R(I/Si). 23 Na is another light nuclei that 
shares most of the characteristics of F albeit with a reduced sensitivity to the SD part, 
see fig. [Ht while 73 Ge is a heavy nuclei that share many of the features of Xe. The ratio 
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Figure 6: Predictions for the ratio of total rates a) Rf/Xc vs Ri/si, b) Rf/Xc vs Ri/n Ge) 
c)RMa/Xe vs Ri/si, d) R[/rsGe vs Ri/si i n different models, same color code as fig. HI 

R F /ri Ge spans roughly the same range than Rf/xc in each model, see fig. Eb- The ratio 
Ri/ 73 Ge just as Ri/Xe does not vary much in either MSSM or LHM, except in the special 
scenarios with much enhanced £ n . In MUED, Rj/Ti Ge goes from 1.5 — 3.5 with high values 
associated with large £/A+, see fig. EH. 

We also considered other combinations of nuclei including those that are primarily sen- 
sitive to spin independent interactions. We found that the rates for R^Ge/Xe-, RAr/Xe, Rat/Cs 
or Rcs/Xe could also give a handle to discriminate models. The results of our scan of the 
parameter space in each of our sample models are displayed in fig. [7] for a representative 
set of pairs of nuclei. As explained before, the rate in n(Xe) increases when £/A + is large 
and the SD contribution important. This increase is driven by so is more important in 
MSSMH than in MUED or LHM. Therefore reductions in the ratio R76Qe/Xe or RAr/Xe 
are larger in MSSMH than in other models, reaching almost a factor 2. On the other hand 
133 Cs is sensitive to £ p thus the rate n(Cs) can be large in MUED models leading to a 
suppression of Rat/Cs while for other models the predictions for Rat/Cs are very similar to 
those for RAr/Xe- Note that in the MSSM, there are two disconnected regions in fig. [7h- 
c. Although the difference is too small to be measured the narrow band corresponds to 
models with a CDM around 50GeV. 

Heavy nucleis can also be used to identify models where tan0 is large. Consider for 
example Rcs/Xe- Both nuclei have similar atomic number and are mainly sensitive to SI 
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Figure 7: Predictions for the ratio of total rates in a) Rp/Xe vs R7e Ge / Xe , b) i?//5, vs 
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interactions, while 133 Cs is also sensitive to £ p and 129 Xe to £„. In models where £ p 3> £ n , 
such as MUED and some of the MSSMQ and LHM scenarios, the ratio Rcs/Xe can be as 
large as 2 while Rcs/Xe ~ 1 in all models where SI interactions dominate and/or £ p w £ n . 

We conclude that in principle with the observation of signals in detectors with different 
materials, including detectors highly sensitive to SD interactions, discrimination of MUED 
from MSSMH and LHM is possible except in a small numbers of scenarios. On the other 
hand the model MSSMQ with light squarks can easily be confused with MUED and LHM. 
We do not attempt to estimate the precision to which the ratio of rates can be measured, 
it is strongly dependent on the specific detector in operation. A more precise analysis 
must await some signals. 

4.4 SI interactions 

The case where only a signal is observed in SI interactions has much less discriminating 
power. Basically the only information that can be used is the total cross-section as well 
as the mass of the CDM particle. We have compared predictions for a^ N for all models 
considered previously including only scenarios where a signal would be seen only in the 
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Figure 8: Predictions for a^ N in scenarios where < 4. x 10 7 pb a) RHNM (pink), 
MUED (red), LHM( green) and b) MSSMH (green) and IDM (black). In b) only models 
that have 0.094 < Vth 2 < 0.136 are included. 

SI channel, that is a^ N > 10~ 10 pb and cr^D < 4 x 10" 7 pb. Here we include also the IDM 
which leads only to a signal in SI interactions. As discussed previously, the predictions for 
the RHNM are always large a^ N > 1.5 x 10~ 7 pb and only a light CDM can be expected, 
see fig. while those of the LHM are much lower 10 _8 pb > a'^ N > 10~ 10 pb. In the 
MUED model, cross sections are very low and the CDM is in the TeV range making it 
very difficult to see a signal. In the IDM model the CDM is either expected to be around 
50GeV with cross sections that span over the full range while for heavier CDM particles 
the predictions do not exceed a few 10~ 9 pb. In the MSSM (only results for MSSMH are 
displayed in the figure, similar results are found for MSSMQ) the whole range of cross 
sections can be expected. For the IDM and MSSM we have imposed the WMAP lower 
and upper limit for the relic density. Allowing other dark matter candidates would lead 
to many more models passing the constraints and cross sections over the whole range for 
any CDM mass. To summarize, with only a signal in the SI channel one could distinguish 
RHNM from LHM, while IDM and MSSMH are often indistinguishable from each other 
and from the previous two models. A DM mass between 500 — 700 GeV is however 
only compatible with the MSSM. Furthermore a heavier DM particle is only compatible 
with MSSM and IDM models while no signal is expected in the MUED model. These 
statements depend crucially on the upper bound that can be set on cr^. For example, if 
no signal is observed at the level of 10 -5 pb we would predict that the SI cross section in 
MUED could reach a^ N ~ 10 -9 pb, an order of magnitude more than what we have used 
in this section. 



5 Conclusion 

We have summarized the predictions for and cr^f in a variety of new physics models. 
We have emphasized the importance of measuring afP to discriminate dark matter mod- 
els although first signals are generally expected in the SI mode. This is because nuclei 
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sensitive to SI interactions provide a measurement of basically one specific combination 
of couplings. On the other hand both amplitudes t; p and £ n can be measured with SD 
sensitive nuclei. Furthermore most particle physics models predict X p « A n , at least those 
where the Higgs is responsible for SI interactions, while £ p /£ n vary over a very wide range, 
from £ p /£ n ~ —1.1 in models dominated by Z exchange to either very small or very large 
values when coloured particles play an important role. To control astrophysical uncer- 
tainties as well as other theoretical uncertainties, we advocate to compare ratios of rates 
measured with different materials. These ratios can be powerful model discriminators es- 
pecially when involving one light nuclei. In particular we have shown that in principle one 
could disentangle the MSSMH and MUED models. Recall that MSSMH is an example of 
a model with a Majorana neutrino where SI/SD interactions are dominated by H/Z ex- 
change while MUED is a model where the same diagrams (quark exchange) contribute to 
both type of processes. With direct detection alone it is much more difficult to distinguish 
MUED models from the LHM scenarios where the new quarks are almost degenerate with 
the CDM as well as with MSSM models with light squarks. This is because in this case 
the SI/SD amplitudes are also both dominated by the exchange of a coloured particle. 
Fortunately if these coloured particles are below the 2 TeV scale, the LHC should be able 
to detect these new particles. 

We have in this analysis concentrated on scenarios with detection rates accessible 
by future CDM detectors. There is however no guarantee that a positive signal will be 
measured for the full parameter space of the models we have considered. This is especially 
an issue for light nuclei. We have also shown explicitly that if the SD elastic scattering 
cross section is below the sensitivity of future detectors, it is much more difficult to identify 
the particle physics model with data from SI direct detection alone, The MSSM, LHM 
and IDM models all predict cr£f in a wide range while a signal should be expected soon 
in the RHNM. In this case no signal is expected in MUED scenarios. A determination of 
the mass of the DM particle will help disentangling some models. 
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